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ABSTRACT 

The newly discovered brown-dwarf eclipsing binary 2MASS J05352184— 0546085 provides a 
unique laboratory for testing the predictions of theoretical models of brown-dwarf formation 
and evolution. The finding that the lower-mass brown dwarf in this system is hotter than its 
higher-mass companion represents a challenge to brown-dwarf evolutionary models, none of which 
predict this behavior. Here we present updated determinations of the basic physical properties 
of 2M0535-05, bolstering the surprising reversal of temperatures with mass in this system. We 
compare these measurements with widely used brown-dwarf evolutionary tracks, and find that 
the temperature reversal can be explained by some models if the components of 2M0535-05 
are mildly non-coeval, possibly consistent with dynamical simulations of brown-dwarf formation. 
Alternatively, a strong magnetic field on the higher-mass brown dwarf might explain its anoma- 
lously low surface temperature, consistent with emerging evidence that convection is suppressed 
in magnetically active, low-mass stars. Finally, we discuss future observational and theoretical 
work needed to further characterize and understand this benchmark system. 

Subject headings: stars: low mass, brown dwarfs — binaries: eclipsing — stars: fundamental parameters- 
stars: formation— stars: individual (2MASS J05352184-0546085) 



Introduction 



Born with masses betwee n the least massive st ars (~ 0.072 M ; Chabrier fc Barafiel 2000 ) and the most 



massive planets (« 0.013 M©; iBurrows et al.ll200lh . brown dwarfs at once extend our understanding of the 
formati on and evolution of both stars and plan ets. In the decade since the discovery of the first brown 



dwarfs (jNakaiima et al 



1995 



Rebolo et al.lll995[ ). the study of star- and planet-formation has increasingly 



trained its attention on these obje cts that, wh ile being neither star nor planet, may provide key insights to 
understanding the origins of both |Basrill2000h . 



Such an understanding must be founded on accurate measurements of fundamental physical properties — 
masses, radii, and luminosities. Unfortunately, the number of objects for which these properties can be 
measured directly is extremely small. Though rare, eclipsing binarie s have long been employed as ideal labo- 
ratories for directly measuring fundamental stellar parameters (e.g. lAndersenlll991l ). The power of eclipsing 
binaries lies in their provision of masses and radii with only the most basic of theoretical assumptions. With 
the addition of atmosphere models the ratio of effective temperatures can also be derived, and luminosities 
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can then be determined directly through the Stefan-Boltzmann law without knowledge of distance. Fi- 
nally, considering the two objects' properties together permits study of the binary as twins at birth whose 
evolutionary histories differ because of their different masses. 

Thus the recent discovery of a brown-dwarf eclipsing binary system in the Orion Nebula Cluster, 2MASS 
J05352184— 0546085 (hereafter 2M0535-05), offers a unique laboratory with whic h to directly and accurately 



test the predictions of theoretical models of brown-dwarf formation and evolution (jStassun. Mathieu. fc Valenti 



20061 , hereafter Paper I). In Paper I, we presented a preliminary analysis of the orbit and of the mutual eclipses 
in 2M0535-05 to directly measure the masses and radii of both components, as well as the ratio of their 
effective temperatures. Our mass measurements reveal both objects in this young binary system to be sub- 
stellar, with masses of Mi = 55 Mj up and Mi = 35 Mj up , accurate to ~ 10%. In addition, from the observed 
eclipse durations and orbital velocities we directly measured the radii of the brown dwarfs to be R\ = 0.67 
R© and R2 = 0.51 R©, accurate to ~ 5%, and representing the first direct measurements of brown-dwarf 
radii. Such large radii are generally consistent with theoretical predictions of young brown dwarfs in the 
earliest stages of gravitational contraction. 

Surprisingly, however, we reported in Paper I that the lower-mass brown dwarf has an effective tem- 
perature that is slightly — but significantly — warmer than its higher-mass companion. Such a reversal of 
temperatures with mass is not predicted by any theoretical model for coeval brown dwarfs. This finding 
has potentially important ramifications for theoretical brown-dwarf evolutionary tracks, and thus for our 
understanding of brown-dwarf formation more generally. 

In this paper, we present supplementary light-curve and radial-velocity measurements (S}2|) which we 
use to refine our determination of the basic physical properties of 2M0535-05 (©. We verify the finding 
of a temperature reversal with mass in this system, and quantitatively compare the empirically determined 
physical parameters with those predicted by theoretical evolutionary tracks. We then briefly explore the 
implications of the temperature reversal in 2M0535-05 for theoretical models of brown-dwarf formation and 
evolution (Sj3J. Specifically, we consider whether the temperature reversal can be explained by dynamical 
brown-dwarf formation scenarios and/or by the physical effects of strong surface fields on young brown 
dwarfs. We summarize our conclusions in Sj5j 



Data and Methods 



2.1. Spectroscopic observations: Radial velocities and spectral types 



We observed 2M0535-05 with the Phoenix spectrograph on Gemini South on eight separate nights from 
December 2002 to January 2003 (first reported in Paper I) and again on 03 March 2 005 (newly r e porte d 



3pon 

here). All nine observations were obtained with the same instrument setup. Following iMazeh et al.l (|2002t ). 



we observed the wavelength range 1.5515-1.5585 fim (central wavelength 1.555 /jm; H band) at a resolving 
power of R ~ 30000 (slit size of 0'.'35). Exposure times were varied between 1.0 and 3.3 hr by the queue 
observers based on sky conditions. The typical exposure time was 2.6 hr (i.e. ~ 0.01 orbital phase; see 
§2.2| , divided into 6 telescope nods in an abbaab pattern for the purposes of sky subtraction and cosmic-ray 
rejection. An obs ervation of the CORAVEL-elodie high-precision radial-velocity standard HD 50778 (K4 III; 
Udrv et alJll999h was obtained immediately before or after each observation of 2M0535-05 to monitor the 
instrument wavelength zero point. In addition, observat ions were obtained of a grid of late-type spectra l 
standards, with spectral types of M0 to M9, selected from lMazeh et al" ( 2002 ') and Mohantv &: Basri ( 2003 ). 
The standard-star observations are summarized in Table [1] 
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The observations were processed using Interactive Data Language procedures developed specifically for 
optimal extraction of Phoenix spectra. The background subtraction step includes logic to compensate for 
variations in the night sky lines due to grating motion or brightness variations between nods. A bootstrap 
procedure is used to determine a wavelength solution from ThNeAr calibration spectra obtained immediately 
before and after the science exposures. First an a priori wavelength dis persion is determin ed by matching 



a Phoenix spectrum of a K giant with the KPNO FTS atlas of Arcturus (jHinkle et al.lll995l) . Next, a small 
wavelength shift (less than 1 km s _1 ) is applied to match the mean wavelengths of the three strongest 
ThNeAr lines as a function of position along the slit. Finally, a multi-Gaussian fit is used to empirically 
determine the wavelengths of the remaining ThNeAr lines. The individual exposures from the different slit 
positions are combined and interpolated onto a uniform dispersion and then continuum normalized using 
sixth-order polynomial fits to the high points in the spectra. Typical signal-to-noise ratios of the extracted 
spectra were ~ 15 per resolution element for the nine observations of 2M0535-05 and ~ 50 per resolution 
element for the standard stars. 

Absolute heliocentric radial velocities of the late- type standards were determined relative to the CORAVEL- 
elodie standard star HD 50778 (see above). To minimize the possibility of introducing systematic velocity 
offsets due to large differences in spectral type between this K4 standard and the late-M standards, we 
applied velocity corrections in a step-wise fashion, determining the velocity shift of each sequentially later 
type standard by cross-correlating it with the standard from the previous step (e.g., MO relative to K4, Ml 
relative to MO, and so on ). The radial velocities reported below are thus formally on the CORAVEL-elodie 



system (|Udrv et alJll999h 



Radial velocities of the tw o components o f 2M0535-05 were determined via the technique of Broadening 
Functions (BFs) described by Rucinski ( 19991) . As discussed in that study, BFs are less prone than simple 



cross-correlation techniques to "peak pulling," in which closely spaced correlation features in double-lined 
binaries can alter the positions of the correlation peak centroids, particularly when the velocity separation 
of the components is on the order of the spectral resolution. The BF analysis requires a radial-velocity 
template spectrum that is well matched in spectral type to the target; we determined that the M6.5 star 
LHS 292 from our grid of standards was the best template for this purpose (see below). 

In seven of the nine Phoenix spectra of 2M0535-05 the BFs show two clear peaks, an example of which 
is shown in Fig. [TJ For these seven ob servations, rad ial velocities and their uncertainties were determined 
using two-Gaussian fits to the BFs (e.g. Rucinski 1999h . We designate the component producing the stronger 



peak in the BFs as the 'primary' and the component producing the weaker peak in the BFs the 'secondary' 
(as we show below, the primary component so defined is also the more massive of the pair). One of the 
remaining two observations was by chance taken during eclipse of the primary (orbital phase 0.75; see Fig. 
[3|) and thus appears as a single-lined spectrum; the BF peak was fit with a single Gaussian and its velocity 
assigned to the secondary. The other observation occurs shortly after eclipse of the secondary (orbital phase 
0.12) and thus appears as a blended spectrum producing a single, broadened BF peak; in this case two 
Gaussians were fit, their initial centroids estimated from the preliminary orbit solution of Paper I and their 
widths fixed based on the widths of the well-separated BF peaks from the other spectra. Radial velocities 
and their uncertainties so measured are summarized in Tableland displayed in Fig.O 

We found that the BF peaks of both components appeared strongest when we used LHS 292, an old M6.5 
dwarf (see Table [lj, as the radial- velocity template in the BF analysis. The BF peaks were ~ 20% weaker 
when we used the M6 or M7 dwarfs as templates, and weaker still when we used earlier or later templates. 
This suggests that the spectral types of both components of 2M0535-05 are also ~M6.5. Indeed, we found 
in Paper I (cf. Fig. 3 in that paper) that the isolated spectrum of the primary component of 2M0535-05 
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very closely matches that of LHS 292, and that the JHK colors of 2M0535-05 imply a spectral type in 
combined light of M6.5±0.5 with negligible reddening. Taken together, the available evidence suggests that 
the components of 2M0535-05 have very similar spectral types of ~M6.5±0.5, corresponding to T r fr « 2700 K 



Slesnick et al. 


2004: 


Golimowski et al. 


2004) 



While our light-curve analysis below ( §2.3[) yields a T c s ratio of ~ 1, reinforcing the conclusion that 
the components of 2M0535-05 are very similar in spectral type, we caution that our assignment of absolute 
spectral types is preliminary and subject to systematic uncertainty. The differences in spe ctral features 



seen i n late-M standard stars are extremely subtle for dwarfs with spectral types M4-M8 (e.g. iBender et al 



20051; IPrato et al.ll2002h . and the low S/N of our 2M0535-05 spectra prohibits a detailed analysis of spectral 



features. Moreover, the low surface gravities and strong magnetic fields of young, low-mass objects can 
bias classification of their spectra, in the sense that t he true spectral types m ay be 1-2 subtypes later than 



inferred from comparison to old, high-gravity dwarfs ijMohantv fc Basrill2003f K Thus, for our purposes here 



we emphasize that the spectral match between 2M0535-05 and the adopted radial- velocity template is suffi- 
ciently good to permit precise radial velocity measurements, and that the spectral types of the components 
of 2M0535-05 are evidently very similar to one another. 

Because the components of 2M0535-05 have such similar spectral types, the relativ e BF peak areas 



reflec t the two components' relative contributions to the total light of the system (e.g. iBavless fc Orosz 



2006). In particular, the primary component apparently contributes approximately 50% more flux than the 
secondary at 1.555 /xm (see Fig.[lJ. This flux ratio provides an important additional constraint for removing 
degeneracies in the determination of physical parameters, as we discuss in §2.31 



The component masses that we determine in H2.3I follow directly from an orbit solution fit to the 
observed radial velocities, and thus the uncertainties in the derived masses (and in other properties that in 
turn depend on those masses) depend sensitively on uncertainties in the radial velocity measurements. The 
mean precisions of the radial-velocity measurements from the BF analysis are formally 1.7 and 1.9 km s" 1 
for the 2M0535-05 primary and secondary, respectively (see <trv column of Table [2]). However, the accuracy 
of these measurements can potentially be degraded by various systematic effects. At the instrument level, 
the stability of the wavelength zero-point is of particular concern. To assess this, we cross-correlated the 
first observation of the CORAVEL-elodie standard star HD 50778 against each subsequent observation of 
this star, providing a measure of the instrument stability over the course of our Phoenix observations. The 
resulting radial velocities of HD 50778 exhibit an r.m.s. scatter of 0.45 km s _1 on the nights that 2M0535-05 
was observed, which is much smaller than the random errors in the individual measurements of 2M0535-05. 
Another possible source of systematic error is spectral-type mismatch between 2M0535-05 and the radial- 
velocity template. However, as discussed above, the radial- velocity template that we used in the BF analysis 
appears to match the spectral types of the 2M0535-05 components very well. Moreover, the components of 
2M0535-05 have very similar spectral types to one another, so that any systematic effects due to spectral 
mismatch with the template should affect both components similarly, in which case only the center-of-mass 
velocity will be affected. 

Finally, the observed radial velocities can be affected by spots on the surfaces of the brown dwarfs that 
cause phase-dependent distortions in the shapes of the spectral lines, and in fact we show in i j2.3.1l that spots 
are clearly present on one or both of the components in 2M0535-05. While the data suggest that this effect 
is small, we have not yet incorporated a physical spot model into our analysis and thus any spot signatures 
remain as potential systematics in the radial velocities. In VI 174 Ori, for example, a young eclipsing binary 
with photometric spot amplitudes similar to those observed in 2M0535-05, we found that the resul ting radial- 



velocity distortions were as large as ~ 1 km s 1 at certain orbital phases (jStassun et al.ll2004at ). Analyses 
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of other spott ed, low-mass eclipsing system s find similar effects; for example, GU Boo sho ws distortions of 
~ 0.5 km s _1 (jLopez-Morales fc Ribasll2005l ) and YY Gem shows distortions of ~ 1 km s _1 (jTorres fc Ribas 
2002). Any spot-induced radial- velocity distortions in 2M0535-05 would need to be 2-3 times larger than 
in these systems t o be comparable to the random measurement errors of ~ 1.5-2 km s _1 (for an example of 



such a system, see lNeuhaeuser et al.lll998l) . Encouragingly the residuals of the radial- velocity measurements 
with respect to the final orbit solution of H2.3\ (xt — 0-7 and 1.2 for the primary and secondary velocities, 
respectively; see Table [2]) do not indicate that systematic effects dominate the radial- velocity measurement 



2.2. Photometric observations: Light curves 

We have continued to photometrically monitor 2M0535-05 intensively with the 1.0-m and 1.3-m tele- 
scopes at CTIO and, as of this writing, possess a total of 2404 J^-band measurements spanning the time 
period 1994 December to 2006 April. The observing campaign to date is summarized in Table [3j and the 
individual measurements are provided in Table 2] We have excluded a small number of observations with 
photometric errors larger than 0.1 mag. Thus, photometric errors on the individual measurements have a 
range of 0.01-0.1 mag, with a mean error of 0.024 mag and a median error of 0.020 mag. 

A period search based on the phase dispersion minimization (PDM) technique of Stellingwerl ( 1978h 



reveals an unambiguous period of P = 9.7795557 ±0.000019 d. The PDM technique is well suited to periodic 
variability that is highly non-sinusoidal in nature, as is the case for most detached eclipsing binaries. This 
updated period is slightly shorter than, but not inconsistent with, the period reported in Paper I (see 
Table EJ). 

In Fig. [3] we show the Ip-band light curve of 2M0535-05 folded on this period. Two distinct eclipses 
are clearly evident and cleanly separated in phase, as is typical of fully detached eclipsing binaries. One 
eclipse — the 'primary eclipse' at orbital phase ~ 0.075 — is notably deeper than the other; this marks the 
time in the orbit when the hotter component is eclipsed which, as discussed in Paper I and below in §3.2| in 
this system is the lower-mass component. 

In addition, the time from primary to secondary eclipse is longer than that from secondary to primary 
eclipse, indicating an eccentric orbit. More quantitatively, the two parameters that determine the shape 
and orientation of the orbit — the eccentricity, e, and the argument of periastron, lu — can be estimated from 
geometrical considerations relating the orbital period, the durations of the ecli pses (6.8 hr and 10.1 hr, 
respectively), and their 6.5-day separation in time (e.g. iKallrath fc Milone1ll999h . from which we estimate 



e « 0.35 and uj ~ 216°. These initial estimates of e and 10 agree well with the more precise values that we 
obtain from detailed light-curve modeling and analysis below ( §2.3| Table [5]). 

In addition to the extensive /p-band light curve measurements, we have also now obtained a set of near- 



infrar ed (JHK) light curves. These are presented and analyzed in a subsequent paper (| Gomez Maaueo Chew et al 



20071 ). but we note here that the same features observed in the Jc-band light curve are also seen at these 



other wavelengths. 
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2.3. Analysis 

As in Paper I, we have performed a simultaneous a nalysis of the Phoen i x rad ial velocities and the J, 



band light curve using the eclipsing-binary algorithms of Wilson &: Devinnevl (jl97ll updated 2005; hereafter 



WD) as implemented in the PHOEBE code of lPrsa fc Zwitterl (|2005l ). WD has become standard in the modeling 



c- 



and analysis of all manner of eclipsing binary systems, having grown in sophistication over the past 35 years 
to include, e.g., phase-dependent projection effects in eccentric orbits, radial- velocity perturbations arising 
from proximity and eclipse effects, apsidal motion, and reflection and limb-darkening effects. In its most 
advanced implementation, the code can also treat the effects of surface spots, model atmospheres, and 
asynchronous rotation, including any attendant gravity brightening and radial- velocity perturbations arising 
from non-sphericity and from phase-dependent variations in the shapes of the components. 

For the present study, we add information from the out-of-eclipse variations in the light curve to model 
asynchronous rotation of the components and as a first step toward accounting for the presence of surface 
spots f £|2.3.ip . As part of our ongoing study of 2M0535-05, we plan to incorporate progressively more ad- 
vanced treatments, including a full investigation of surface spots and inclusion of state-of-the-art brown-dwarf 
model atmospheres. Here, as in Paper I, we use simple blackbody spectra in our light curve models. While 
brown-dwarf spectra are known to exhibit strong departures from simple blackbodies, we show below that 
this effect is unlikely to significantly alter our principal findings; in particular, the reversal of temperatures 
with mass that we find in 2M0535-05 is probably not the result of improperly modeled spectra. On the 
other hand, understanding the temperature reversal may ultimately require a proper treatment of surface 
spots and their influence on the time- and aspect-dependent emergent flux from the system. 



2.3.1. Spot-related variations in the light curve 



To begin, we ran an initial WD fit to the light-curve and radial- velocity data (Tables [2] and [4j using as 
initial values the system's orbital and physical parameters from Paper I, together with the updated orbital 
period from this study f §2.2|) . The best-fit parameters and their formal uncertainties are summarized in 
Table [SJ The reduced x 2 of this fit is large (xt — 3.10), and the r.m.s. of the residuals relative to this fit 
(0.033 mag) is accordingly larger than expected given the typical error on the photometric measurements 
fi )2.2p . This suggests an additional source of variability in the light curve — spots, for example — not included 
in the model. 

Surface spots are most commonly manifested as low-amp litude, periodic, roughly sinusoidal variations 
superposed on the eclipse light curve (e.g. lStassun et al 1120044 - We searched for such l ow-amplitude , periodic 
variations in the light curve of 2M0535-05 by applying a Lomb-Scargle periodogram |Scarglelll982h analysis 
to those portions of the light-curve data obtained outside of eclipse. We performed this periodogram analysis 
separately on each epoch of light -curve data (se e Tabic [3l), as spot- related variations often show evolution in 
amplitude and phase with time ( Stassun et al. 1999t Torres fc Ribasl 20021 Stassun et al. 2004a ). 



We find statistically significant periodic signals in all but the first two epochs, these earlier epochs 
possessing relatively few measurements. While the amplitudes of these periodic signals vary slightly from 
epoch to epoch, they are always small (Am < 0.06 mag, peak-to-peak). Importantly, the periods of these 
signals are consistent across epochs, with a value of P S pot ~ 3.3 d. Spot-related variability observed in other 
eclipsing binaries is similarly ch aracterized b y photometric amplitudes of a few percent and periods that are 
constant over many epochs (e.g. lRibas1l2003h . most likely reflecting the rotation period of one or both binary 
components. The implication here is that at least one of the brown dwarfs in 2M0535-05 is nearly always 
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spotted, and rotates with P rot w 3.3 d. 



With a rotation period that is shorter than the orbital period by a factor of (almost exactly) 3, the 
rotational and orbital angular velocities of 2M0535-05 have evidently not yet become pseudo-synchronized 
with one a nother (wh ere the rotational angular velocity is synchronized to the orbital angular velocity at 
periastron; iHutl Il98lh . which for the eccentricity of 2M0535-05 would give P ro t/-Forb ~ 1/2- Instead, the 
rotation of at least one component in 2M0535-05 is super-synchronous. Such super-synchronicity is probably 
not surprising given the extreme youth of 2M0535-05 (see below) and that young brown dwarfs tend to be 
rapid rotators (Moha ntv &: Basri 2003); tidal effects have likely not yet had sufficient time to synchronize 
the system (e.g. iMeibom fc Mathieull2005h . 



A more detailed discussion of rotation in 2M0535-05 is beyond the scope of this paper. For our present 
purposes, we proceed to adopt the observed P ro t = 3.3 d in our final WD fit as a constraint on the oblateness 
of both 2M0535-05 components. In addition, we rectify the /p-band light curve by subtracting a sinusoid of 
P = 3.3 d, with the amplitude and phase determined separately for each epoch of data, as described above. 
The aim of this rectification procedure is to "remove" the spot signal from the light curve and to thereby 
allow the wd model to achieve goodness of fit without introducing additional spot-fitting parameters that 
are poorly constrained at present. We caution, however, that this procedure is not equivalent to modeling 
the spots physically. The inclusio n of additional light curves at multiple wavelengths into our analyses 



(jGomez Maaueo Chew et al.l 120071 ) will ultimately be required to constrain the physical properties of the 
spots and to permit full modeling of their effects on the observed light curves and radial velocities. Bearing 
these caveats in mind, we proceed with our analysis below using the rectified ip-band light curve, and revisit 
the possible influence of spots in §3.21 



2.3.2. Determination of physical parameters 

The final wd solution was determined by simultaneously fitting the Phoenix radial velocities and the 
Jc-band light curve, rectified as described above. The resulting orbital and physical parameters of 2M0535- 
05 are summarized in Table El the fit to the radial velocities shown in Fig. [21 the fit to the light curve 
displayed in Fig. [31 and the system geometry illustrated in Fig. [U With the light curve rectified, the WD fit 
is now excellent; the r.m.s. of the residuals is 0.021 mag, comparable to the typical error on the photometric 
measurements f N2.2p , and the goodness-of-fit is correspondingly very good, xl — 1-035. 

A comparison of the orbital and physical properties of 2M0535-05 derived from the rectified and unrec- 
tified light curves (Tabled reveals mostly minor, statistically insignificant, differences in the fit parameters. 
However, a few parameters do show differences that are larger than their formal uncertainties. For example, 
both e and ui differ by ~ 2a (although these parameters are strongly correlated with one another, so their 
uncertainties are not independent), suggesting that the rectification procedure is not entirely free of small, 
systematic couplings to some parameters. Thus, while a physical treatment of spots in the WD model will 
ultimately be required to achieve full accuracy in the determination of some parmaters, the effects of spots 
on the most fundamental physical parameters that we derive for 2M0535-05 — the orbital parameters, in par- 
ticular, and thus the component masses and radii — are probably insignificant. Moreover, if other parameters 
remain susceptible to spot-modeling effects, these effects are evidently very subtle. 

We note that there exists an alternative WD solution to the one presented in Tabled with nearly identical 
best-fit parameters except that the component radii are reversed. Formally, the goodness-of-fit of this alter- 
native solution is inferior, with xt = 1-062, but it is nonetheless acceptable. This particular type of degener- 
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acy i s a general feature of eclipsing-binary solutions (see, e.g., IStassun et al.ll2004at lLopez-Morales &: Ribas 
2005). Since only the component radii are significantly different in the two solutions, the degeneracy can 



be broken by adding the additional constraint of the luminosity ratio, which is generally very different in 
the two solutions. In this case, the preferred solution from Table [5] gives an ii-band luminosity ratio of 
L1/L2 = 1.63 ± 0.13, whereas the alternative solution gives L1/L2 = 0.69 ± 0.06. The luminosity ratio 
predicted by the alternative solution — in which the secondary is both warmer and larger, and thus more 
luminous, than the primary — is clearly inconsistent with the ratio inferred from the observed spectra, which 
imply L1/L2 ~ 1.5 in the H band (see Fig. [IJ. 



3. Results 

3.1. Masses and Radii 

With masses of Mi = 0.0569 ± 0.0046 M Q and M 2 = 0.0358 ± 0.0028 M , the components of 2M0535- 
05 are proven to be bona fide brown dwarfs. Moreover, with radii of R\ — 0.674 ± 0.023 R© and R2 = 
0.485 ± 0.018 R©, these brown dwarfs are more comparable in their physical dimensions to young, low-mass 
stars than to old brown dwarfs, consistent with their being very young (see below). Importantly, these mass 
and radius measurements are distance independent. Moreover, systematic errors do not appear to dominate 
over random measurement uncertainties (~i j2.ll) . Thus the masses and radii of 2M0535-05 reported here are 
accurate to 8% and 3%, respectively. In M3.4I we make use of these measurements for an initial examination 
of the predictions of theoretical brown-dwarf evolutionary models. 



3.2. Reversal of temperatures 

The physical properties that we have determined for the brown dwarfs in 2M0535-05 are surprising in 
at least one important respect: The less- massive secondary is hotter than the primary. Specifically, we find 
T e ff,2/T e ff,i = 1.064 ± 0.004 (Tabled, which follows from the relative depths of the eclipses and bolometric 
corrections from the model atmospheres (blackbodies in this case), with small corrections for differences 
in occulted areas that occur due to the mildly eccentric orbit (Fig. [4|) . This finding is highly statistically 
significant. Such a reversal of temperatures with mass is not predicted by any theoretical model for coeval 
brown dwarfs — in which temperature increases monotonically with mass — and thus merits closer scrutiny. 
Here we consider the possible effects of non-blackbody atmospheres and surface spots. 

Brown dwarf spectra deviate significantly from an ideal blackbody, due primarily to the strong wave- 
length dependence of molecular opacity. The corresponding redistribution of flux into or out of any specific 
wavelength depends on temperature, potentially affecting the ratio of eclipse depths. However, the primary 
and secondary of 2M0535-05 have very similar effective temperatures f fl2.1[) . so both components should 
have similar deviations from a blackbody. 

To check this quantitatively, we have examined up-to-date synthetic models of brown-dwarf spectra 
for temperatures appropriate to the comp o nents of 2M0535-05. In particular, we use the solar-metallicity 
CONE0 atmosphere models of Allard et al. ( 2001 ). As the masses and radii of 2M0535-05 imply log g « 3.5 



1 At the temperatures of the components of 2M0 535-05, the DUSTY atmo sphere models of these authors are very similar and 
would work equally well for our purposes. See also lMohantv fc Basril l|2003l l. 
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for both components, we use the COND models with this logg value also. The COND models incorporate 
recent opacities, including ~ 500 million molecular lines, as well as dust formation and condensation at very 
low temperatures. (For a more detailed discussion of these models, and their applicabi l ity to young brown 
dwarfs in particular, we point the reader to the extensive discussion in lMohantv fe Basril (|2003) ) . The results 
of our analysis are summarized in Fig. [5] We find that, while dramatic departures from simple blackbodies 
are clearly present in the model spectra, the ratio of surface brightnesses in the I band is consistent with 
that predicted from simple blackbodies to within 2%, which implies a T c g ratio that is consistent with the 
blackbody value to within 0.5%. 

As an additional check, we calculated the T e ff's of the components of 2M0535-05 independently at 
1.25, 1.65, and 2.2 /im. We use the ratio of surface fluxes at 0.8 /im (measured directly from the ratio of 
eclipse depths in the Ic-band light curve; §2.3|) and the radii (measured directly from the eclipse durations 
and the orbital velocities; Tabled]), together with an ass umed distance to the Or ion Nebula Cluster of 450 
pc, empirical bolometric corrections fr om the literature ( Golimowski et al.ll2004f ). and the observed JHK 
magnitudes at these wavelengths from ICarpenter et al.l |200ll ): K = 13.58 ± 0.02, J - H = 0.640 ± 0.015, 
and H-K = 0.385 ± 0.015. We derive temperatures of 7\ = 2730 K and T 2 = 2875 K (1.25 fxm), Ti = 2675 
K and T 2 = 2860 K (1.65 fim), and T x = 2680 K and T 2 = 2795 K (2.2 /im). In all cases, the derived T off 's 
are consistent with a spectral type of ~M6.5, and in all cases the reversal of temperatures persists. 

Another potential contributor to non-blackbody spectra is the presence of hot or cold spots (akin to solar 
plage or sunspots) on the brown dwarfs, leading to spectra that derive fro m a combination of temp eratures. 
Cool spots have now been found to be present on many brown dwarfs (e.g. Scholz fc Eisloffel 2005) and may 
thus be expected on one or both brown dwarfs in 2M0535-05 as well. Indeed, we find a clear periodic signal 
in the out-of-eclipse portions of the /c-band light curve, with a period of P S pot ~ 3.3 d and a semi-amplitude 
of ~ 0.03 mag f ^2.3.1|) ; almost certainly this is a rotationally modulated spot signal. 

Such spot signals are commonly observed in the light curves of close eclipsing binaries, though modeling 
these effects is often a largely cosmetic exercise with little effect o n the resulting stellar properties. For 
example, in our analysis of the young eclipsing binary VI 174 Ori (jStassun et al. I l2004ah . we successfully 
reproduced the observed out-of-eclipse light-curve variations by including spots in the light-curve model, but 
the derived system parameters were altered only negligibly as a result (see also lTorres fc Ribasll2002t iRibas 
20031 : lLopez-Morales fc Ribasll2005h . 



On the other hand, accurate modeling of spot effects can be crucially important for the proper inter- 
pretation of certain systems. W UMa stars of the so-called "W type" are a particularly good case in point. 
A defining characteristic of these systems is that the deeper eclipse corresponds to the occultation of the 
physically smaller an d lower-mass sec ondary, similar to what is observed in 2M0535-05. One interpretation 
advanced early on bv iRucinskil (|1974l ) is that the secondary is hotter than the primary (typically by ~ 5%), 
perhaps due to thermal effects arising from mass - and energy-transfer in these contact systems. However, 
an alternate interpretation (e.g. Eaton et al. 1980) is that cool spots on the primary, if preferentially located 
along the the eclipsed latitude and more-or-less uniformly distributed in longitude, can have the effect of 
lowering the surface brightness of the eclipsed regions on the primary during transit by the secondary, thus 
making the eclipse of the primary shallower. Such a model is generally very delicate in the details, and 
requires careful arrangement of the putative spots in order that they produce only mild (< 0.05 mag) out-of- 
eclipse variations i n the light cu rve; the inferred spot configurations can resemble a "leopard print" pattern 
(e.g. lLinnel]jll99lh . 



m some cases 



We are now experimenting with more sophisticated light-curve models that include the effects of spots. 
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In the meantime, we emphasize here that the ratio of eclipsed surface brightnesses implied by our current 
light-curve model is robust; the higher-mass primary really does radiate less per unit area than the secondary, 
at least at those surface elements that are eclipsed by the secondary. Thus, in what follows, we continue to 
explore the implications of a temperature reversal in 2M0535-05. 



3.3. Physical association with the Orion Nebula star-forming region: Evidence for youth 



Strong evidence for the physical association of 2M0535-05 with the Orion Nebula Cluster, and hence of 
its youth, is provided by its center-of-mass velocity. The observed value of 7 = 24.1 ± 0.4 km s _1 (Table [5]) 
is within 1 km s _1 of the systemic radial velocity (25 ± 1.5 km s -1 ) of kinematic members of this active 
star-forming region (jstassun et al.lll999t IsTcilia-Aguilar et alj|2005h . 



In addition, we can derive a distance to 2M0535-05 by comparing the total system luminosity to its 
observed flux. To calculate the luminosity, we use the directly measured radii (Table [5]) together with the ef- 
fective temperatures and apply the Stefan-Boltzmann relation, L = 4ttR 2 <7bT^ s . The adopted spectral type 
of M6.5±0.5 for the 2M0535-0 5 primary (32T1) implies T eSA = 2715 ± 100 K based on recent calibr ations of 
the T G ff scale for brown dwarfs (jMohantv k Basrill2003l ; ISlesnick et alj|2004t iGolimowski et al.ll20041) , though 



as mentioned earlier, systematic uncertainties in the T e g scale may be as large as ~ 200 K. From the measured 
temperature ratio of T e s, 2 /T eS ,i = 1.064 ± 0.004 (Table HJ, we obtain T cff , 2 = 2820 ± 105 K. This then gives 
component luminosities of L\ = 0.0223 ± 0.0034 L Q and L2 = 0.0148 ± 0.0023 L , and a to tal system lumi- 
nosit y of L = 0.0372 ± 0.0060 L©. Adopting an apparent magnitud e of ttik = 13.58 ± 0.02 (ICarpenter et al 



20011 ) and bolometric corrections appropriate for the observed T e ff 's (jGolimowski et al.ll2004f ) yields a derived 
distance to 2M0535-05 of 456 ± 34 pc, assuming no extinction. This distanc e determination is consistent 
with the distance to the Orion Nebula of 480 ± 80 pc (jGenzel k Stutzkil [l989h . 



The Orion Nebula Cluster is extremely young, with an age that has been estimated to be just lj^ 
Myr ( Palla k Stahler 1999 : Hillenbrand! 1997 ). With a center-of-mass velocity and distance that are both 
consistent with membership in this cluster, 2M0535-05 is probably also very young, likely having formed 
within the past few Myr. In addition, as discussed in Paper I, the JHK colors of 2M0535-05 place an upper 
limit on the extinction of Ay < 0.75, and thus limit the amount of remnant material available to the brown 
dwarfs for ongoing accretion. Given the eclipsing nature of the system, any disk material would necessarily 
be seen edge-on and would thus produce a large amount of extinction and reddening. Thus, while the colors 
of 2M0535-05 are formally consistent with a small amount of interstellar extinction and/or a small amount 
of remnant disk material, the currently observed masses are unlikely to change significantly over time. These 
brown dwarfs will likely forever remain brown dwarfs. 



3.4. Comparison to theoretical brown-dwarf evolutionary tracks 



While we expect that our ongoing analysis of 2M0535-05 will further improve the accuracy of the 
measured system parameters, the present accuracy is sufficient to permit an initial examination of theoretical 
brown-dwarf e volut ionary models. Here we consider two o f the m ore widely used sets of tracks: those of 
Baraffe et all (jl998 .fl and those of b'Antona k Mazzitell3 |l997i updated 1998). These models differ in 



2 We use the tracks with convection mixing-length parameter a = 1.0, as recent analyses of young, low-mass stars and 
brown dwarfs suggest that such low-efficiency convection best matches the observed physical properties of these objects (e.g. 
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their treatment of brown-dwarf atmospheres and of energy transport (convection) in brown-dwarf interiors, 
as well as in the choice of initial condit ions (for a more in-depth discussion of differences in these models, 
see 



Siess et all 12000: iBaraffe et alj |2002) 



Fig. [6] compares the observed radii, effective temperatures, and luminosities of 2M0535-05 with the 
values theoretically predicted by these models for brown dwarfs with masses of Mi — 0.0569 ± 0.0046 M Q 
and M 2 = 0.0358 ± 0.0028 M Q (Table EJ). 

We note first that, generally speaking, the agreement between the observed and theoretically predicted 
properties of 2M0535-05 is quite good. The models predict that these very young brown dwarfs should, at an 
age of ~ 1 Myr, be significantly larger, warmer, and more luminous than their older counterparts — and that 
is in fact what we see. Indeed, between 1 and 100 Myr, a brown dwarf with a mass of 0.057 M© (the mass 
of the 2M0535-05 primary) is predicted to shrink by 500%, cool by several hundred K (or, equivalently, go 
from an M spectral class to an L spectral class), and dim by 1.5 orders of magnitude (Fig. [S]). The recently 
measured radius (R = 0.12 R^) o f the old and low-mass (M = 0.09 M©) star in the eclipsing binary system 
OGLE-TR-122 ijPont et alJbood) confirms that stars with near-brown-dwarf masses have very small radii 
(~ 1 Rjup) when they are old. Thus, the fact that 2M0535-05 comprises young brown dwarfs that are both 
large and luminous — and even simply that they are of M spectral class, as predicted — is a testament to the 
generally good predictive power of current theoretical models of brown dwarfs. 



At a more detailed level, the models of iD'Antona &: Mazzitellil (|l997l ) predict radii and luminosities 
that are mor e cons istent with the observed values. The difference is most pronounced in the radii; the 
Baraffe et al.l (|1998l ) models under-predict both radii by ~ 10%. This finding of under-predicted model 
radii is similar in sense and magnitude to that found from recent efforts to der ive the physical proper ties of 
low-mass stars and brown dwarfs through detailed modeling of their spectra (jMohantv et al.ll2004al lbr). In 
addition, while both sets of models p erform reasonably well with re spect to the observed luminosities, the 
agreement is somewhat better for the ID'Antona fc Mazzitellil (|1997l ) model s; it appears possibl e that more 



(|1998I ) models for 



accurate measurements will reveal an under-prediction of luminosity by the IBaraffe et al 
the lower-mass secondary. 

But perhaps more importantly, the reversal of temperatures in 2M0535-05 remains puzzling; the rela- 
tionship between effective temperature and mass is predicted by both sets of models to be monotonic for 
brown dwarfs of the same age. However, this expectation disappears for two brown dwarfs o f differ ing age. 
The observed effective temperatures can be seen as consistent with the ID'Antona fc Mazzitellil (|1997l ) models 
if the primary is taken as being modestly younger than the secondary (At 0.5 Myr). Indeed, the tem- 
peratures, radii, and luminosities of 2M0535-05 all remain in marginally good agreement with these models 
for such an age difference, at a mean age of 1 Myr (Fig. |7|), if we also adjust the observed T e s scale cooler 
by 70 K. Such a shift is well within the current systematic uncertainty of the brown-dwarf T e g scale (£ !2.ip . 
Thus, a question raised by our findings is whether current brown-dwarf formation theory can accommodate 
a scenario in which two brown dwarfs — that are part of the same, very young, binary system — can have 
sufficiently different ages to allow for reversed temperatures as we have observed in 2M0535-05. 



IStassun et aUl2004al; iMathieu et~al"ll2006l ; iTorres et aTbood) . 
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4. Discussion 

4.1. 2M0535 05 as a case study in dynamical brown-dwarf formation? 



In truth, we have very little un derstanding; about how bina ries with components as close as those 



in 2M053 5-05 are formed (see , e.g., bonnelll EoOlL iTohlinel Eool . (2M0535-05 is, after PP1 15 in the 



Pleiades ( Basri fc Martini 1999), the shortes t-period brown-dwarf binary system yet discovered). With _ 
few exceptions (e.g. iTohline fc Durisen 2001 ). fission seems to be ruled out. In situ mechanisms involving 
dynamic cloud fragmentation and disk fragmentation have not yet proven successful in creating such close 
binaries, but remain to be fully developed. In particular, the role of orbital migration in the presence of a 
massive circumbinary disk needs to be considered. Certainly, if the two brown dwarfs formed together as a 
close binary, we have essentially no a priori expectations for whether they should be observed to be coeval 
to within 0.5 Myr. 

Alternatively, recen t theoretical work ( Reipurth fc Clarke] 2001 ) , as well as detailed numerical simula- 
tions ( Bate et al. 2002al b: Bate fc Bonnelll 2005h . suggest that dynamical interactions may be integral to the 
formation of brown dwarfs. The argument is essentially that strong gravitational interactions in multiple- 
body encounters provide a feasible mechanism for disrupting the accretion process, and thereby preventing 
the accumul ation of mass by ob jects that would otherwise have become stars. This hypothesis remains under 
debate (e.g. iMaxted fc Jeffries! 120051 ). Still, it is tempting to speculate that the components of 2M0535-05 
did not form together as a binary, but rather formed separately — with the primary forming later — and then 
were later married through a dynamical interaction. In such a scenario, it is possible that the resulting binary 
system — comprising two objects that were not originally formed together — may exhibit seemingly peculiar 
characteristics, such as the observed reversal of effective temperatures, that in fact reveal the non-coeval 
nature of the system. 

There are at least two specific scenarios involving multiple-body interactions that might pertain to the 
origin of 2M0535-05. One involves a low-mass binary pair that interacts with a more massive third body from 
elsewhere in the cluster. Simulations show that, if the binary is not simply broken apart by the encounter. 



the lo wer-mass member of the binary is ejected and replaced by the massive incoming object (|Bate et al 



2002a|). The resulting binary would then likely consist of non-coeval members. A serious concern with this 
scenario in the context of 2M0535-05 is that it is unlikely in three-body encounters within clusters that the 
most massive object would have a mass of only 0.06 M©. 

Thus, a second scenario is that several objects formed in the fragmentation of a small molecular core, 
forming an unstable multiple system. The consequen t rapid dynamical evolution of the system may have 
both terminated accretion and formed a hard binary |Bate et al.ll2002bh . In such a fragmentation scenario 
the relative core-collapse times and evolutionary zero points might easily vary by 0.5 Myr, of order the age 
difference needed to explain the temperature reversal of 2M0535-05 in the context of some evolutionary 
models (Fig. [7]). 

The position of 2M0535-05 might be taken as further evidence that dynamics have been important in its 
history. With a projected separation of 2.8 pc from the center of the Orio n Nebula, 2M0535-05 is s i tuate d 
more than 10 core radii from the center of the Orion Nebula Cluster ( Hillenbrand fc Hartmann 1998); 
perhaps the pair was ejected from the cluster cent er during a dynamical encounter. The one-dimensional 



velocity dispersion in the cluster is ~ 2 km s 1 ( van Altena et al. 1988t Hillenbrand fc Hartmannl 1998t 



Stassun et al. 1999 : Sicilia-Aguilar et al. 200ot) . so to reach its current position in 1 Myr, if ejected from the 
cluster core, 2M0535-05 would need a somewhat higher-than-usual velocity of ~ 3 km s _1 in the plane of 
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the sky. However, the measured center-of-mass radial velocity of 2M0535-05 does not deviate significantly 
from the cluster velocity. An alternative interpretation for the position of 2M0535-05 at the outskirts of 
the cluster may be that it is not a member of the Cluster proper, but rather a me mber of the more widel y 
distributed population of young stars in the region surrounding the Orion Nebula ( Warren fc Hesser 1978j) . 
The binary's proper motion is the critical measurement needed to establish whether 2M0535-05 was ejected 
from the cluster center. 



4.2. Missing physics in theoretical brown-dwarf evolutionary tracks? 



An alternative explanation is that the evolutionary tracks are deficient with respect to some critical 
physical ingredient (s). For example, the presence of strong magnetic fields on one or both brown dwarfs in 
2M0535-05 could be affecting energy transport and thereby altering their physical structure. 

I ndeed, recent analy ses of several young, low-mass eclipsing binaries fe.g. lStassun et al. L j ICovino et al 



2004 iTorres et al.ll2006l ) indicate systematic discrepancies in the models. In particular, the observed effective 
temperatures are cooler than expected, and the observed radii larger than expected. These discrepancies 
are especially pronounced among the most magnetically active stars (as traced by X-ray emission and other 
proxies). One possible interpretation is that strong magnetic fields inhibit energy transport in these other- 
wise fully convective stars, resulting in a decr ease of the surface temperatur e and a corresponding increase 
in radius so as to radiate the same total flux jMontalban fc D'Antorialbood) . Such an interpretation would 
be consistent with, and would help explain, the emerging observational e vidence for suppre ssed convection 
in young, low-mass stars (e.g. Stassun et al. 2004a ; Mathieu et al. 20061 Torres et al. 20061 and references 
therein). 

In this context, the temperature reversal in 2M0535-05 might be taken as evidence for a strong magnetic 
field on the higher-mass primary that is causing a sufficient decrease in its surface temperature to make it 
effectively cooler than the lower-mass secondary. By inference, the secondary would be interpreted as being 
less magnetically active. The observational eviden ce is strong for magnetic activity in young, low-mass stars 
and b rown dwarfs of early- and mid-M type (e.g. iMohantv et al.ll2002t IStassun et al.ll2004i iPreibisch et al. 
20051). Moreover, the evidence in fact show s a marked decline in magnetic activity at very late M types 



(e.g. Gizis et al. 2000 ; Mohantv et al. 20021 ) . suggesting that brown dwarfs with roughly the mass of the 
2M0535-05 secondary and below are not capable of generating strong fields. The idea of a magnetically 
active primary would also be consistent with the primary being heavily spotted (see £ !3.2j) . 



If we are to explain the anomalously low effective temperature of the 2M0535-05 primary in this way, 
we should then also exp ect its radius to be larger tha n theoretically predicted, as the one effect goes hand 
in hand with the other |Montalban fc D'Antonalbood ). As discussed above ( 33~4l). whether the 2M0 535-05 
radii agree with theoretical predictions depends on one's choice of model. The Baraffe et al. ( 19981 ) tracks 
do indeed suggest over-sized radii in 2M0535-05 (for both c omponents), though this may simply reflect the 
truncation of those tracks at 1 Myr. On the other hand, the lD'Antona fc Mazzitellil (|1997l) tracks agree with 
the observed radii very well (Fig. [6]), with no need to invoke missing physics. 
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5. Summary and Conclusions 

2M0535-05 is the first known eclipsing binary system comprising two brown dwarfs. Satisfying both 
kinematic and distance requirements for physical association with the young (~ 1 Myr) Orion Nebula Cluster, 
2M0535-05 provides the only direct, accurate measurements of the fundamental physical properties of newly 
formed sub-stellar objects. The masses that we measure are accurate to ~ 10%, the radii accurate to ~ 5%, 
the ratio of effective temperatures accurate to ~ 1%, and all are distance independent. As such, 2M0535-05 
represents an important benchmark for theoretical models of brown-dwarf formation and evolution. 

Encouragingly, we find that current brown-dwarf evolutionary tracks are, broadly speaking, successful 
in predicting the fundamental physical properties of these young brown dwa rfs. More quantitatively , of the 



two sets of theoretical models considered here, we find that the models of iD'Antona fc Mazzitellil (| 19971 ) 
yield mass-rad ius and mass- l umino sity relationships that best agree with the empirically determined ones. 
The models of Baraffe et al. ( 19981) predict radii and luminosities that are 1.5-2(7 smaller than the observed 



values. 

However, the reversal of component effective temperatures with mass in 2M0535-05 is unexpected and 
unexplained. We have considered here two possible interpretations of this intriguing result. The first is that 
the components of 2M0535-05 are m ildly non-coeval, with the high er-mass primary being ~ 0.5 Myr younger 



than the secondary The models of ID'Antona fc Mazzitellil (|1997| ) are in fact consistent with the observed 



temperature reversal for such an age difference. A second hypothesis is that strong magnetic activity on the 
primary is inhibiting convection, and thereby lowering its surface temperature. 

Neither of these interpretations is wholly satisfying, and neither is obviously discreditable. Binary 
formation theory is largely silent on the subject of coevality at the level of ~ 0.5 Myr, and theorists have 
long warned the star- formation community about the limited applicability of evolutionary track chronometry 
at such early ages (the model zero-points being arbitrary in most cases). In addition, while the observational 
evidence is strong that young brown dwarfs can be magnetically active, the effects of magnetic fields on 
brown-dwarf structure and evolution have yet to be consistently modeled or fully understood. 
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Table 1. Standard stars 



Name 


SpTy 


Exp. time (s) 


Ref. 


GJ 328 


MO 


180 


Mazeh et al. (2002) 


GJ 382 


Ml. 5 


180 


Mazeh et al. (2002) 


GJ 447 


M4 


120 


Mohantv & Basri (2003) 


GJ 406 


M6 


120 


Mohantv & Basri (2003) 


LHS 292 


M6.5 


3600 


Mazeh et al. (2002) 


LHS 3003 


M7 


900 


Mohantv & Basri (2003) 


GJ 3655 


M8 


3600 


Mohantv & Basri (2003) 


BRIB 1507-0229 


M9 


3600 


Mohantv & Basri (2003) 
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Table 2. Radial velocity measurements of 2M0535-05 



HJD a 


Phase b 


R.V. C 


CRV 


(O 


-C) d 






km s _1 


km s _1 


km 


s" 1 






P 


rimary 








2452623.74805 


0.75 










2452624.75701 


0.85 


9.59 


1.41 




1.34 


2452625.72555 


0.95 


2.00 


2.22 




1.32 


2452626.65154 


0.04 


9.84 


1.75 




-2.79 


2452649.72938 


0.40 


38.02 


1.41 




1.30 


2452650.67281 


0.50 


34.21 


1.14 




0.31 


2452655.74512 


0.02 


6.61 


2.56 




•0.98 


2452656.70824 


0.12 


25.58 


1.28 




-0.84 


2453432.57920 


0.45 


37.09 


2.11 




1.62 










= 0.7 




Secondary 


2452623.74805 


0.75 


29.69 


1.17 




2.50 


2452624.75701 


0.85 


48.53 


1.73 




•0.41 


2452625.72555 


0.95 


61.35 


2.07 




0.40 


2452626.65154 


0.04 


42.43 


2.97 




0.43 


2452649.72938 


0.40 


4.18 


1.24 




0.37 


2452650.67281 


0.50 


5.07 


2.00 




3.21 


2452655.74512 


0.02 


50.04 


2.35 




0.05 


2452656.70824 


0.12 


22.58 


1.31 




2.45 


2453432.57920 


0.45 


6.64 


2.64 




0.85 



Xl = 1.2 



a Heliocentric Julian Date. 

b Orbital phase, relative to ephemeris of Table [5] 

c Heliocentric radial velocity. 

d Residual relative to WD solution (see text). 
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Table 3. Summary of Ic-band time-series photometry 



Obs. 


Telescope 


HJD 


range a 




iV obs b 


1 


KPNO 0.9-m 


49699. 


.70 


- 49703 


,93 


40 


2 


WISE 1.0-m 


49698. 


,35 


- 49714 


,42 


47 


3 


WIYN 0.9-m 


52227. 


.78 


- 52238 


.01 


102 


4 


WIYN 0.9-m 


52595, 


,75 


- 52624 


.95 


144 


5 


SMARTS 0.9-m 


52622 


,57 


- 52631 


.51 


122 


6 


SMARTS 1.3-m 


52922. 


,73 


- 53081 


.57 


347 


7 


SMARTS 0.9-m 


53011 


,57 


- 53024 


.77 


205 


8 


SMARTS 1.3-m 


53280. 


,74 


- 53340 


.73 


230 


9 


SMARTS 1.0-m 


53373. 


,56 


- 53386 


.79 


184 


10 


SMARTS 1.3-m 


53403. 


,53 


- 53445 


.60 


338 


11 


SMARTS 1.3-m 


53646 


.82 


- 53727 


.69 


204 


12 


SMARTS 1.0-m 


53719. 


.59 


- 53727 


,83 


190 


13 


SMARTS 1.3-m 


53745. 


.64 


- 53846 


.48 


251 



a Range of Heliocentric Julian Dates (2400000+) . 
b Numbcr of observations. 
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Table 4. Differential 7 c -band light curve of 2M0535-05 



HJD a 


Am b 


Cm 


Obs. c 


2453686.70525 


0.006 


0.020 


11 


2453686.71293 


0.016 


0.054 


11 


2453687.76841 


-0.004 


0.052 


11 


2453687.77617 


-0.006 


0.052 


11 


2453687.78389 


0.020 


0.029 


11 


2453687.79158 


-0.002 


0.029 


11 


2453688.73523 


-0.001 


0.020 


11 


2453688.74295 


-0.008 


0.061 


11 


2453688.75066 


0.022 


0.020 


11 


2453688.75835 


0.030 


0.020 


11 


2453689.77018 


0.262 


0.025 


11 


2453689.77783 


0.220 


0.020 


11 


2453689.79309 


0.237 


0.062 


11 


2453690.71340 


0.008 


0.020 


11 



a Heliocentric Julian Date 

b Differential Ic magnitude (arbitrary 
zero-point). 

c Source of measurement (see Tabled]). 



Note. — The full table is available in 
the electronic version of the Journal. A 
portion is shown here for guidance regard- 
ing its form and content. 



Table 5. Orbital and physical parameters of 2M0535-05 





Paper I 


This 

Unrectificd a 


study 

Rectified" 1 


Orbital period, ± [clj 


q 77Q(^oi -i- 0000/19 


9.779556 ±0.000019 b 


1 Ul JJLT Ido II UI1 I -DL-oot-IIIdll jt,dl I. 1Q 


9001 Sfi^fi^O -I- OOOOQ^. 


2001.863903 ±0.000160 


2001.863765 ±0.000071 




n qooc. _l n no^o 


0.3354 ±0.0049 


0.3276 ±0.0033 


Orientation of periastron, uj [°j 


91 E /I _|_ 1 1 
Z10.4 zc 1.1 


219.2 ± 1.4 


217.0 ±0.9 


Semi-major axis, asinz [a.u.] 


o^cis -i- n 001 
u.uoyo zt u.uuiu 


0.0406 ±0.0016 


0.0406 ±0.0010 


Center-of-mass velocity, 'y [km s ^] 


o/i 1 _i_ n a 


24.1 ±0.4 


24.1 ±0.4 


1\/Tpcg T*nHr_ n — l\/f n I A/f 
IVIdoD IdllU, y — JVJ2 / iW^ 


n ^oe _i_ n 01 8 


0.622 ± 0.022 


0.631 ±0.015 


Tntal rvift^Q ( A/f . -1- A/f«^ Qin^ 7 [M^l 

lULcll Illdoo, lil-i^ 1 J*J2/' 3 ^^^ ^ 


n n&&0 4- 007fi 

U.UOOU Z1Z U.UU 1 u 


0.0932 ±0.0111 


0.0932 ± 0.0073 


Tn i n £i 1~ i pin 7 [°1 

_LIH_.IlIld LIUII , t I 


00 _i_ n 9 

OO.O —1— U . __j 


89.4 ±0.3 


89.2 ±0.2 


_r iniidi y odiij-diiijjiii uul, i \ \ Kin o 




18.37 ± 1.01 


18.49 ±0.67 


Secondary semi-amplitude, K.% [km s _1 ] 




29.55 ± 1.24 


29.30 ±0.81 


Primary mass, M\ [M©] 


0.0541 ± 0.0046 


0.0575 ±0.0069 


0.0572 ± 0.0045 


Secondary mass, M2 [M©] 


0.0340 ± 0.0027 


0.0358 ±0.0043 


0.0360 ± 0.0028 


Primary radius, R\ [R©] 


0.669 ±0.034 


0.673 ±0.037 


0.675 ±0.023 


Secondary radius, R2 [R©] 


0.511 ±0.026 


0.485 ±0.029 


0.486 ±0.018 


Primary gravity, log gi 




3.62 ±0.14 


3.62 ±0.10 


Secondary gravity, log #2 




3.54 ±0.14 


3.54 ±0.09 


Effective temperature ratio, TijT\ 


1.054 ±0.006 


1.062 ±0.006 


1.064 ±0.004 



'WD solutions based on fits to the unrectified and rectified /c-band light curve (see §2.31) . 



b Uncertainty in the period is from a phase dispersion minimization ( Stellingwera 1 1 978 ) analysis on the Jc-band light 
curve (see £|2 . 2[) : the period is held fixed in the WD fit and its uncertainty propagated into the uncertainties of derived 
quantities. 
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Fig. 1. — Determination of the radial velocitie s of the primar y and secondary components of 2M0535-05 
using the Broadening Function (BF) analysis of lRucinski ( 1999f) . The example BF shown here (solid line) is 



from the spectrum of H JD 2452649 (see Table [5]) , where the primary and secondary are at near maximum 
velocity separation. An M6.5 dwarf was used as the radial- velocity template (see M2.1[> . The dashed line shows 
a two-gaussian fit to the BF, from which velocity ccntroids and their formal uncertainties are determined. 
The ratio of the peak areas from this fit is 1.6, indicating that the primary dominates the light of the system, 
contributing ~ 60% more flux than the secondary at the wavelength of these observations (1.555 /im). 

Fig. 2. — Radial velocity measurements of 2M0535-05. The individual radial velocity measurements from 
Table [2] are plotted (primary measurements in green, secondary measurements in red), folded on the orbital 
period and phased to the time of periastron (see Table [5|) . The solid lines are WD models based on a 
simultaneous fit to the radial velocity measurements and the rectified Ip-band light curve (see Table |4] and 
Fig. [3]). Distortions in the model curves near phases 0.075 and 0.75 are due to the brief occultations of 
the approaching and receding limbs of each component when it is eclipsed. Residuals are shown at bottom. 
Note: This figure appears in color in the electronic version of the journal. 

Fig. 3. — /c-band light curve of 2M0535-05. The individual photometric measurements from Table |4] are 
plotted, rectified using sinusoidal fits to the out-of-eclipse portions of the light curve (see Ej2.3[) . folded on 
the orbital period and phased to the time of periastron (see Table [5]) . The solid line is a WD model based 
on a simultaneous fit to the rectified Jc-band light curve and the Phoenix radial velocities (see Table [2]) . 
Residuals are shown at bottom. Insets show detail around primary and secondary eclipses, which in this 
system correspond to the eclipses of the secondary and primary components, respectively. The r.m.s. residual 
is 0.02 mag, comparable to the mean photometric error. The reduced % 2 of the fit is 1.035. Note: this figure 
appears in color in the electronic version of the journal. 

Fig. 4. — The geometry of the orbital plane of 2M0535-05 in the rest frame of the primary (more massive) 
brown dwarf is illustrated, to scale, at primary eclipse (orbital phase 0.075; see Fig. [3]). The brown dwarfs are 
represented by filled circles, their radii also to scale (the adopted rotation period of P ro t = 3.3 d constrains 
the oblateness of both components to be less than 0.5%). The position of periastron (orbital phase 0.0) 
is indicated by a star symbol, and small arrows indicate the direction of the secondary's orbit about the 
primary. The observer is to the left, as indicated by the arrow. Note that, at this orbital phase corresponding 
to the deeper eclipse, the secondary (less massive) brown dwarf is the one eclipsed. 



Fig. 5. — We use the COND atmosphere models of lAllard et al.l (|200lh to check whether the non-blackbody 



spectra of brown dwarfs may explain the unexpected reversal of effective temperatures with mass in 2M0535- 
05. Upper panel: COND model spectra at the effective temperatures of the primary (green) and secondary 
(red) components of 2M0535-05 are compared with simple blackbodies of the same temperatures (dashed 
lines). Dotted vertical lines demarcate the approximate bandpass of the Iq filter used for our light-curve 
observations. Bottom left: The flux ratio of the COND models in the upper panel is shown relative to the 
blackbody flux ratio as a function of wavelength in the Iq band. Also shown is the ip-band filter bandpass 
profile (dotted curve) and the bandpass- weighted mean COND-to-blackbody flux ratio (dashed line), which 
is 1.02 in this case; that is, the COND flux ratio is within 2% of the blackbody value. Bottom right: Same as 
bottom left, except showing the T e ff ratio implied by the COND flux ratio. The bandpass-weighted value is 
1.069 (dashed line), within 0.5% of the value of 1.064 ± 0.004 found in our WD analysis f£ )2.3p using simple 
blackbodies. Note: This figure appears in color in the electronic version of the journal. 



Fig. 6.- 



Comparisons of observations with theoretical models of young brown dwarfs. The measured radii, 
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effective temperatures, and luminosities of 2M0535-05 (symbols with error bars) are compa red to the values 
predicted by two sets of theoretical models ( Baraffe et al.ll 1998L D'Antona fc Mazzitelli|[l997f) of young brown 



dwarfs. Solid curves show the predicted evolution from 0.1 to 100 Myr for brown dwarfs with masses equal to 
those measure d for the prim a ry (gr een) and secondary (red) brown dwarfs in 2M0535-05. (The theoretical 
calculations of Baraffe et al. ( 19981 ) do not extend to ages less than 1 Myr.) Dashed curves bracketing the 
solid curves represent the la measurement uncertaintie s in those masses. Obs e rved and predicted values 
are generally in good agreement, particularly with the iD'Antona fc Mazzitellil (|1997l ) models. Note that 
the effective temperature of the primary is p redicted by both sets of mode ls to be warmer than that of the 
secondary at any particular age, but that the lD'Antona fe Mazzitellil (jl997l ) models allow the primary to be 
cooler than the secondary if it is sufficiently younger (see also Fig. [7]). Note: This figure appears in color in 
the electronic version of the journal. 



Fig. 7. — Same as Fig. [6j except with the 2M0535-05 components arbitrarily separated in age by 0.5 Myr 
(ti = 0.75 Myr, r 2 = 1.25 Myr). This demonstr ates that the temperature rev ersal observed in 2M0535- 
05 may in fact be marginally consistent with the ID'Antona fc Mazzitellil (|1997l) models if the 2M0535-05 
primary is taken to be slightly younger than the secondary. Note: This figure appears in color in the 
electronic version of the journal. 
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